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We show that the derivative of the logarithm of the average characteristic polynomial of a diffusing Wishart 
matrix obeys an exact partial differential equation valid for an arbitrary value of N, the size of the matrix. 
In the large N limit, this equation generalizes the simple Burgers equation that has been obtained earlier for 
Hermitian or unitary matrices. The solution through the method of characteristics presents singularities that we 
relate to the precursors of shock formation in fluid dynamical equations. The 1 IN corrections may be viewed as 
viscous corrections, with the role of the viscosity being played by the inverse of the doubled dimension of the 
matrix. These corrections are studied through a scaling analysis in the vicinity of the shocks, and one recovers 
in a simple way the universal Bessel oscillations (so-called hard edge singularities) familiar in random matrix 
theory. 
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I. INTRODUCTION 

Recently, two of us have argued [1, 2] that Dyson's orig- 
inal idea H to interpret the distribution of eigenvalues of a 
random matrix as the equilibrium limit of a random walk per- 
formed independently by each of the matrix elements, can be 
exploited in order to get a new insight on some properties of 
the spectrum of eigenvalues. We have considered the two sim- 
plest, albeit non-trivial examples: additive and multiplicative 
matrix-valued Brownian motions. In both cases, the evolving 
spectra aregoverned by complex, non-linear diffusion equa- 
tions |H Ulil HI • The solutions of these equations, based 
on the method of complex characteristics, exhibit singularities 
that we may associate to the precursors of shocks (which we 
refer to as 'pre-shocks' jg]). Interestingly, for finite values of 
N, corresponding to a non-zero "spectral viscosity" v, = 
these shocks are accompanied by universal, oscillatory phe- 
nomena. In the case of additive Gaussian random walk, the 
oscillations are those of an Airy function. In the case of the 
multiplicative matrix valued-diffusion yj], the corresponding 
evolving matrices are unitary, and an additional phenomenon 
occurs when, due to the compactness of the support of eigen- 
values (unitary circle), the left and right Airy shocks merge 
into a more complicated structure described by a Pearcey 
function. This universal dynamics has been observed, for ex- 
ample, in the behavior of Wilson loops in simulations of Yang- 
Mills theory with a large number of colors (8|]. In this case, 
the role of the time is played by the area of the Wilson loop 
and the size of the matrix is determined by the number of col- 
ors. The critical oscillations of the preshock are associated to 
the onset of non-perturbative effects of the Wilson loop, i.e. 
they signal the transition to the strong coupling regime. 
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In this paper, we further develop this fluid dynamical per- 
spective for the family of random matrices known as the 
Wishart ensemble [9]. Originally formulated as the multivari- 
ate analogue of the gamma distribution, the Wishart ensemble, 
plays, till today, a dominant role in multivariate statistical data 
analysis, including such modern fields as image processing, 
genetic chips, financial engineering, wireless telecommunica- 
tion and several others 11 1 CHI . Complex Wishart ensembles are 
also closely related to chiral random matrix models, defining 
the universal regime of spectral properties of Dirac operators 
in Euclidean non-abelian gauge theories fiUl . 

We start by formulating the random walk of complex 
Wishart ensembles in terms of a Smoluchowski-Fokker- 
Planck equation for the probability distribution of eigenval- 
ues. From this equation, after proper rescaling of time, and us- 
ing standard techniques of statistical mechanics, we recover in 
the large N limit, a complex non-linear equation of the Burg- 
ers type for the resolvent (whose imaginary part yields the av- 
erage spectral density). We solve this equation by the method 
of complex characteristics, and recover the well-known spec- 
tral properties of the ensemble (Bronk-Marcenko-Pastur equa- 
tion) lfl2jfl3ll . We also make a link to the Voiculescu Hill free 
random variables theory. We proceed then to the analysis of 
finite N effects in order to capture the universal behavior that 
is anticipated in the vicinity of the singular points associated 
with the edges of the spectrum and that correspond to pre- 
shocks in the Burgers equation. In order to do so, we proceed 
as in the prequels of this paper Q2|, |4f] , and focus on the time- 
dependent average characteristic polynomial, more precisely 
the derivative of its logarithm, which, in the large N limit sat- 
isfies the same Burgers equation as the resolvent. An exact 
nonlinear partial differential equation can be obtained for this 
object, that is valid for arbitrary N, and that therefore allows 
us to study the oscillatory behavior in vicinity of the shocks 
that is visible for large but finite N. In particular, at one of 
the edges, we obtain the Bessel equation, whose solution de- 
fines the universal behavior of the spectral functions near the 
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origin of the spectrum, i.e. so-called "hard edge" singularity. 
We conclude the paper with the study of the evolution of in- 
verse characteristic polynomials. The Cauchy transforms of 
the associated monic orthogonal polynomials fulfill the same 
nonlinear Burgers-like partial differential equations. 



II. THE RANDOM WALK OF WIS HART MATRICES 

We consider a N x N random matrix of the form L = 
K'\t)K( f) belonging to the complex (J3 — 2) Wishart ensemble. 
The matrix K is rectangular MxN (M > N) and with complex 
entries. We assume that the elements of K evolve in 'time' ac- 
cording to dKij(t) = bf)(t) + ibf)(i), where b { })(i),bf)(i) are 
real numbers defined by two independent Brownian walks: 



with 



and 



(^f)4fV)) = crW<5(f-f'). 



(i) 



(2) 



(3) 



This is sometimes referred to as a Laguerre process, an ana- 
logue of the non-colliding squared Bessel process IU5I1 . The 
real valued equivalent (bf^\t) — 0, /? = 1) is called the Wishart 
process, and was first studied in lfl^[l7ll . 

Following Dyson y|], we shall study the behavior of the 
eigenvalues of L as a function of time. Note however that 
since we are not interested in the equilibrium distribution, we 
do not include any restoring force on the motion of the eigen- 
values. That is, the dynamics is entirely given by the ran- 
dom walks described above. The random matrix K admits a 
singular value decomposition K = UkV, with U, V unitary, 
and k = diag(/ci , • • • , kjv) rectangular. The eigenvalues of L, 
which we call At, are related to the non-singular values of K 
by A, = /c?. Furthermore, it is worth mentioning, that the pairs 
-Ki and k, are the nonzero eigenvalues of the complex chiral 
matrix: 



W = 





K 



Note, moreover, that the M x M ('anti- Wishart') matrix L a = 
KK^ has the same set of non vanishing eigenvalues as the 
NxN matrix L. In addition, both L a and W possess M-N = v 
zero eigenvalues. The easiest (and physically most transpar- 
ent) way to obtain the time evolution of the k, 's, is to use sec- 
ond order perturbation theory [3]. A simple generalization of 
the calculation presented in Ref. 11911 yields: 



j8(v+l)-l , P 



K*e> 



1 1 

H 



2/Cj 

(SkiSkj) = S,j 8t. 
To get the evolution of the Ai's, note that 

A! i = (ki + 5Ki) 2 = Aj + iK-fiKi + 5k}, 



8t, 



(4) 
(5) 



so that, upon averaging, (with 5A t = A' t - Aj) 

(SAi) = 2Ki{5Ki) + (6k]). (6) 
Using this relation and Eq. (O one then easily obtains 



(6Ai)=l3 



St, 



{6 Aid A j) =AA i 6 ij dt. 



(7) 



Eqs. (0]i and (0 can be also derived via Ito calculus |2j 
and written in the following stochastic differential forms: 



dKi = bi(f) + ~ 



v+1- 



2 



1 



1 



Kj K j K( H~ Kj 



and 



dA i = 2^ i b i (f)+J3 



v + 1 + 2A t ^ 



jm Ai ~ A J 



>dt, (8) 



dr. (9) 



The latter can be interpreted as an equation for the evolution 
of N, non-intersecting, aka "vicious", walkers confined to the 
positive part of the real axis. It has a physical realization as, 
for example, a model for nonintersecting fluctuating interfaces 
in thermal equilibrium |[2lll . Recently, similar diffusion pro- 
cesses for the Wishart ensembles were also generalized for 
arbitrary value of parameter j3 € (0, 2l j22ll . 

All the statistical properties of the eigenvalues can in princi- 
ple be derived from the probability P(A\,A2, ■■■ , A^, t) to find 
the eigenvalues near the values A\, ■ ■ ■ ,Ay at time t. By us- 
ing Eq. (0 or Eq. (0, we may write a Smoluchowski-Fokker- 
Planck equation for P. It takes the following form: 



m Ai - A ' 



(10) 



In this paper, we shall consider the specific case where the 
random walk ([]]) is initiated at the origin, i.e., K,j(0) = for 
all f s and fs. In this case, it is easy to verify that the solution 
of Eq. GOll is 



P^Nt-^YliAj-AifYlA 
where TV is a normalization constant. 



!(v+i)-i 



(11) 



EQ. DYNAMICS OF THE RESOLVENT IN THE LARGE N 
LIMIT 

In this section, we analyze the bulk properties of the spec- 
trum of eigenvalues by deriving an equation for the average 
density from Eq. (fTOb . This is possible in the limit of large 
matrices. Since we are dealing with rectangular matrices, care 
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must be taken in taking this limit. Let us define the rectangu- 
larity r = N/M and recall that we already set v = M - N. The 
limit M — > oo, N — > oo can be taken either at r fixed, in which 
case v — > oo, or at v fixed, in which case r — * 1. Most of the 
time we shall work at fixed r. 



A. The partial differential equation for the Greens function of 
the Wishart random martix 

We start by recalling some definitions. The averaged den- 
sity of eigenvalues is defined by: 

/N N 
Y\dA k P(.A u ---,A f dY i s ( A - A d = 
k=\ 



(12) 



One defines similarly the 'two-particle' density 

p(A,p,t) = (YYj 6iA - Al)6 (»- A ->i) ■ (13) 
\ 1=1 fat) ' 

These densities are normalized as follows 

j dAp(A, t) = N, j dAdpp{A,p, t) = N(N - 1). (14) 
Furthermore, observe that 

with J denoting the principal value of the integral. 

The equation for p (A, t) is easily obtained by multiplying 
Eq. ( [Tol l by Yn=\ 5 (A - A[) and integrating over all N eigen- 
values. One gets 

dp(A,t) d 2 r 1 [ / 1\ 

J^=2- 2[ Ap(A, t)]+ p\N(l--)-l 

d d \ Cp(A,Li,t) 



(16) 



This equation simplifies in the limit of very large matrices. 
To see that, let us set p (A, p) = p(A)p(p) + p con (A, p), where 
p con (A, p) is the connected part of the two-point density, ex- 
pected to be of order 1 /N as compared to the factorized con- 
tribution p (A) p (p). Next, we simultaneously rescale the time, 
so that r = Nt, and change the normalization of the densities: 



p(A)=Np(A), p(A,p)=N(N-l)p(A,p). 
This results in the following equation 



(17) 



dp(A,r) (1 \ dp(A,T) 
dr \r dA 



<A 



yy 'J A-p \ N 



d 2 p(A,T) 
OA 1 



4-pdp(A,r) 



N 



dA 



H dA J A-p H 



(18) 



where on the right side of the equal sign, we have gathered 
terms that disappear in the limit N,M — > oo with r fixed. 

At this point it is convenient to introduce the usual resolvent 
(sometimes also called Green's function) 



dp 



Pip, t) 



(19) 



z-p 

The imaginary part of G(z - A- ie, r) yields the average spec- 
tral density, whereas the real part is the Hilbert transform of p 
(<H [p (A)] = J j^dp). In the limit N,M —> oo, with r fixed, 
we then obtain the following closed equation for p (A, r): 

d P (A,T) / \\ dp(A,T) 

3t P \ r) 8A + 

d 

-2/3A— \p (A, T)<H\p (A, T )]} - 2p P (A, T)9i\p (A, r)] . 



(20) 



One can take its Hilbert transform to obtain: 

m\p(A,T)\ I _ i\ &h\pu,t)] 

8t P \ rj dA 



8A 

+f3[p(A,T)] 2 -/3{<H[p(A,T)]} 2 . 



dA 



(21) 

where we have used that for f(x) = \p(xfH [p(jc)]}, 
1ri\xf{x)\ = xfH [f(x)]. By combining the two equations 
above, we derive the equation for the resolvent that is the ana- 
logue of the Burgers equation obtained earlier J2, Ht] for the 
GUE and CUE ensembles: 

d T G(z, t) = j8(l- --2 Z G(z, r))d z G(z, t)-/3G 2 (z, t). (22) 
r 

This equation has been obtained in a slightly different form 
in Ref. [23], using other techniques. To establish the corre- 
spondence with the equation derived in Ref. 12311 . we need to 
rescale t — > Then Eq. (l22l transforms into 

d T G(z,r) = -d z G(z,r) 

+r (d z G(z, t) - 2zG(z, f)d z G(z, r) - G 2 (z, r)) , (23) 



in full agreement with 12311 . This nonlinear partial differential 
equation, governing the evolution of the resolvent of a freely 
diffusing Wishart matrix, will be solved in the next subsection 
using the method of (complex) characteristics l23ll . as we did 
in our earlier works (U, 01 • 



B. Solution with complex characteristics 



It is convenient to rewrite Eq. d23l obeyed by G as follows, 



dG 8G , 
(1 - r + 2rzG) — + — + rG 2 = 0. 
dz dr 



(24) 



from which the equations for the characteristic lines (parame- 
terized by s) are easily obtained 



dz 

— = 1 - r + 2rzG, 
ds 



(25) 
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as 



as 



(26) 



(27) 



We choose to solve these equations with the same initial con- 
dition that leads to Eq. ( fTTT ). namely G (z, r = 0) = 1/z. This 
translates into z(s = 0) = zo,r(i = 0) = 0, G(s = 0) = 1/zo- 
Solving the last two equations gives then respectively s = t 
and 



G = 



1 



rr + zo 



We are therefore left with: 
dr 

whose solution reads: 



r + 



2rz 
rr + ZQ 



z = ^1 + — j (zo + rr) . 

One can then use this result to eliminate zo in Eq. 
yields the following implicit equation for G(z, t): 



1 



z = 



G(z,t) 1-rrGfcT) 



(28) 



(29) 



(30) 



which 



(31) 



The construction just outlined supposes that the mapping 
between z and zo is one-to-one, that is, it can be inverted. This 
is the case except at points where dz/dzo = where a sin- 
gularity occurs, that may be interpreted as the formation of a 
"preshock" in the Burgers equation. The singular points are 
determined by 



dz 1 r n 

— = 1- -jT = 0, 

dzo zi 



zqc -±yfh 



(32) 



corresponding to z c = (1 ± -\[P) 2 t. These singularities oc- 
cur precisely at the edges of the spectrum, as we shall verify 
shortly. Note that the left singularity (for positive z) originates 
from characteristics starting at negative zo. 
The solution of Eq. PTT ) reads explicitly 



(r-l)r + z- Vz 2 -2zr(l+r)+T 2 (T^7)2 

G(z,r) = , (33) 

2rrz 

where the (minus) sign in front of the square root is cho- 
sen so that the resulting spectrum density is positive. G(z, t) 
matches its initial condition when r — > 0. We have also that 
G(z, t) behaves as 1/z at large z and fixed r, as appropri- 
ate since the spectrum is bounded. By taking the imaginary 
part of G(z - A - ie, t), one recovers the well known Bronk- 
Marcenko-Pastur formula for the level density: 



y/(A-C- T )(c + T-A) 

2nArr 



(34) 



The spectrum is localized in the interval [c_r, c + r] with c ± = 
(1 ± -\fr) 2 , as anticipated from the study of the preshocks. 



C. The partial differential equations for related ensembles 

A similar treatment can be applied for the ensembles 
closely related to the Wishart ensemble. For instance, the re- 
solvent for the chiral matrix is defined as 



(w, t) = 



1 



N + M 



Tr- 



1 



W(T) 



(35) 



It is connected to its Wishart equivalent by the following trans- 
formation (see e.g. O): 



G(z) = 



r- 1 r+ 1 



2rw 2 2rw 



'(w), 



(36) 



where w = z. By inserting this change of variable into 
Eq. d23l ), one gets 



1 +r 



l+r 



d T g + gd w g 



0. 



(37) 



For M — N, i.e., r = 1, this equation reduces to the same 
Burgers equation that was obtained for the diffusing Hermitian 
matrix in ||2|] . 

We can additionally relate the Green's function of the 
Wishart ensemble to that representing random anti- Wishart 
matrices. We have OHIM]: 



1 / 1 



1 



+ rG(z), 



(38) 



a relation which also holds in the case when the matrices 
evolve in time (the term 1/z corresponds to the zero modes 
of the matrix, which are of 'kinematical' origin, reflecting the 
rank of the matrix). The analog of Eq. ( 1231 for G a is 

d T G a {z,T) = {\-r-2zG a {z,T))d z G a {z,T)-G 2 a {z,r). (39) 

The Green's function for the chiral matrix can be easily ob- 
tained from d33l by using the relation (f36t between G(z, r) and 
g(w, r). One gets 



g (w, t) = 



w 2 - -^/(w 2 -c_t) (w 2 -c + t) 



(r + 1)tw 

Similarly, the anti- Wishart Greens function reads 

(1 -r)r+z- Vz 2 -2zr(l+r) + T 2 (l-r) 2 



G a (z,T) 



2tz 



(40) 



(41) 



Note the pole at z = 0, which reflects the presence of the zero 
modes that we have already mentioned. The residue, which 
reads \—r= ^jf-, accounts exactly for the ratio of the number 
of zero modes to the total number of eigenvalues. 



D. Connection to the Voiculescu approach 

It is interesting to observe that Eq. OTI ) that emerged from 
analyzing the complex characteristics of Eq. d23l can also 
be obtained from the free random variable formalism lfl4ll . 
The cornerstone of this approach is the so-called R-transform, 
which plays the role of the logarithm of the characteristic 
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function in the classical probability calculus. The R-transform 
is additive under free convolution, alike the logarithm of the 
characteristic function, generating additive cumulants for the 
convolution of independent probability distributions. The R- 
transform is related to the Green's function by 



R[G(z)] + l/G(z) 



(42) 



or equivalently G[R(z) + 1/z] = Z. Thus, modulo the shift 
1 lz, the R-transform is the functional inverse of the Green's 
function. 

The R-transform for a static Wishart random matrix is 
known (see e.g. [25]). It reads R st (z) = 1/(1 - rz). In order 
to get the R-transform of the time-dependent Wishart matrix 
L(r), we note that the time evolution of the probability dis- 
tribution P(t) is equivalent to a linear rescaling with time of 
the eigenvalues in a stationary probability distribution func- 
tion P st (see e.g. Eq. (fTTTi). We may therefore write 



Gfcr) = i(Tr — l — 
N\ z- L(t) 



PIT) 



1 / 1 

-(Tr 

N\ z 



= -G st (-)- 
tL s , I Psi t t 



(43) 



Then, by using the definition of the R-transform, Eq. ( 142b . 
and the explicit form of the R-tranform of the Wishart matrix 
recalled above, we get R(z, r) = t/? s ,(tz). We see therefore 
that the quadratic equation OTb emerging as the solution of 
the complex characteristics equation is simply Eq. ( l42b for the 
time-dilated R-transform for the Wishart ensemble. 



IV. EXACT EQUATION FOR THE AVERAGE 
CHARACTERISTIC POLYNOMIAL AND ITS SCALING 
LIMITS 

As it has been already emphasized in 0], it does not ap- 
pear possible to push the simple, and physically intuitive, de- 
scription of the level density given in the previous section, 
beyond the large N limit, since the 1 /N corrections are diffi- 
cult to calculate. In order to study these finite size corrections, 
we shall then proceed as in J2l, and use the technique of or- 
thogonal polynomials. After a suitable generalization of the 
relevant polynomials to the time-dependent problem, we ob- 
tain for these polynomials exact equations. In particular the 
logarithm of the average characteristic polynomial is found to 
obey an exact equation, which in the large N limit coincides 
with the equation for the resolvent that was discussed in the 
previous section. This allows us to study universal behavior 
in the vicinities of the edges of the (large N) spectrum, albeit 
not of the spectral density itself but of the logarithm of the 
average characteristic polynomial. 



A. Time dependent characteristic polynomials 

For the chosen initial conditions, the average characteris- 
tic polynomial associated with the diffusing Wishart matrix is 



equal to a certain time dependent monic orthogonal polyno- 
mial El 



<det[z-L(r)]> = M Ar (z,T). 



(44) 



In the static case the relevant orthogonal polynomials are the 
generalized Laguerre polynomials [28]. These are defined as 



(45) 



and satisfy the orthogonality relation 



' e- x ^L a n {x)L a m {x) = S nm r( " + " +1) . (46) 



The following recursion relation 

nL°{x) = (or + 1 - x)l£\ (x) - xL a n + _j(x), (47) 
and the differentiation property 



— k L°{x) = (-!)*£££(*), 



(48) 



will also be useful. 

For the dynamical case, we need to construct time depen- 
dent, monic, polynomials that are orthogonal with respect to 
the measure defined by the solution of the Smoluchowski- 
Fokker-Planck equation (ITTb for /3 — 2. It is easy to check 
that the following polynomials satisfy these requirements 

M a n {x,r) = (-T) n n\L a n {^, (49) 
with the orthogonality condition given by 

T^) ' M a n {x,T)M a m {x,T) 

= 5 nm (t) 2 " +1 T(n + a + V)n\ . (50) 

The analogues of Eqs. (l47T i and ( 1481 read respectively 

nM a n (x, r) = n[x- r(a + 1)] M"*{ (x, r) + 

-n(n-l )txM°^1 (x, t) , (51) 

and 



Jo 



—rM a „(x,T)= 

dx k (n — k)\ 



M^(X,T). 



(52) 



In order to derive the differential equation satisfied by the 
polynomials M" (x, t) we first observe that 

%-M a n {x, t) = -M a n (x, t) - ™«jc, t). (53) 

OT T T OX 

Then, using ( Bil l, we eliminate the first term on the right hand 
side of this equation, which yields: 

d T M a n (x, t) = -xd x Ml(x, t) - (1 + a) d x M a n {x, t). (54) 

At this point, we rescale the time in order to be consistent with 
our earlier calculations, namely we set: r — > (with f3 = 2). 
We also set a — v. Finally, since we are interested mostly in 
the average characteristic polynomial, we focus on n — N. We 
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thus obtain the following exact partial differential equation for 
the characteristic polynomial: 

d T Ml{z, t)=~ [zd zz Mliz, r)+(l +v) 8M N {z, r)\ . (55) 

It is easy to show that this equation yields Eq. (|23]» for the 
resolvent in the large N limit. To do so, let us recall first that 
d. In <det [z - L(t)]) = NG(z, r) in the large N limit. We then 
define 



f N (z,T)=^d z ln[M v N (z,T)] 



(56) 



(the inverse Cole-Hopf transform), and obtain from Eq. (IB3T l 
the following equation for 



d T f N + r(2z,f N dJ N +f%) + (1 - r)d z f N 
= ~{2dJ N +zd z J N ). 



(57) 



In the large N limit, the right hand side can be dropped, and 
one recovers Eq. d23l for fy(z, t) = G(z, t). 

The exact equation $55[ that we have obtained will be used 
in the following subsection in order to analyze the 1 /N cor- 
rections of the characteristic polynomial in the vicinity of the 
edges of the spectrum. 



B. Characteristic polynomial at the edge of the spectrum 

We are interested in the asymptotic scaling of the averaged 
characteristic polynomial at the edge of the spectrum, where 
the average eigenvalue spacing is of order N~ s . The power 
6 is determined by the singular behavior of the spectrum (in 
the large N limit) near its edge. It is easily verified that a 
power law behavior ~\z — z c \ a yields 5-1 /(a +1). We shall 
then write z = Z c (j) + N~ 6 s, and study the behavior of the 
characteristic polynomial, or its inverse Cole-Hopf transform, 
as a function of s in the large N limit. We shall set 



fy (z«r(T) + N~ s s, t) * A(t) + AT^(s, r), 



(58) 



where the function %(s, t) remains finite in the large N limit, 
and y — 1 - 5 (for a singularity ~ \z- z c \ a , J — a I (a + 1)). 

Two distinct limiting behaviors arise: the first one, which 
occurs at the right shock for any r and at the left one for r + 1 , 
is referred to as 'soft edge scaling' ; the second, obtained at the 
left shock for r = 1, i.e. at the origin, is known as 'hard edge 
scaling'. We shall examine successively these two cases. 



1. The soft edge 

We know the solution / M (z, t) of Eq. d57b in the limit 
N,M—>oo and fixed r. This is given by Eq. d33l ). i.e., 



,, , (r-Yyr + z- y/(z - z L )(z - z R ) 

U(z,t) = , (59) 

2rrz 

where zl = t(1— -\fr) 2 and zr — t(1 + -\/r) 2 are the positions 
of the left and right edges of the spectrum, respectively. As 



announced above, we either probe the close vicinity of the left 
z = Zl + N~ s s (r + 1) or the right z — Zr + N~ s s traveling 
shocks. At those points the average eigenvalue spacing is pro- 
portional to N~ 2 ^ and so 6 = I, corresponding to a square 



root singularity (a = 1 /2). Our ansatz 
following form: 



has therefore the 



tib 2 (r- 1)t + z . 1 
ft fc/* +N~s,t) = — 2» +at**(j,t), (60) 

Inserting it into d57b and keeping the dominant terms as Af — > 
oo, one obtains: 



+ 2xd s x + d ss x = 0, 

z* 



(61) 



9 1 9 1 

where we define z* — {-z^r^^z^r^-T} depending on the edge 
we look at. This equation can be written 



d s 



X + dsX ~ - 



o. 



which is easily integrated to yield 



X +dsX~-+ g(f) = 0, 



(62) 



(63) 



where g(j) is an arbitrary function of r. In terms of <p(s, t) de- 
fined as ^(i, t) = d s ln 4>(s,t) (inverse Cole-Hopf transform) 
one obtains 



(i,T) + |g(T)--|0(i,T) = O. 



(64) 



The shift of variable s = y+g(r)z*, and the redefinition if/(y) = 
(p(y + g(r)Zm), transform this equation into the equation for the 
Airy function 



y 

d yy 4t{y,T) i/<y,T) = 0, 



The solution reads: 



<p(s, r) = Ai 



s - g(T)z, 



(65) 



(66) 



In order to find the unknown function g(r) we match the large 
s asymptotics of d60b with the Green's function (133b at \z\ — > 
z L/R ■ This way (see the Appendix) we find that g(r) = 0. 



2. The hard edge 

The spectrum of the infinite size L(r) matrix, for any t > 0, 
"touches" the origin only if the rectangularity r - 1. This is 
therefore when we can observe the hard edge universal scal- 
ing. To this end we rewrite d57b in the following form: 

d r f N +r(2zf N d z f N +fl) = ~ [{2 + v)d z f N +zd zz f N } . (67) 

In the limit of M and Af going to infinity with v fixed and finite 
(so that r — > 1), we have: 



/oo(Z,T) 



z - ^z 1 - 4tz 
2tz 



(68) 
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In the vicinity of the origin, the singularity is of the inverse 
square root type, so that the average eigenvalue spacing is pro- 
portional to N~ 2 , or 5 = 2. We set then z = A^s and Eq. d58l l 
takes the following form: 

MN~ 2 s,t) = ±+Nx(s,t). (69) 

Inserting d6*9l into (l67l > and taking the large N limit with r — > 1 
and v fixed, one obtains: 

sd ssX + 2s X d sX + (2 + v) 8 sX + X 2 = 0. (70) 
We now have 

d s [sd sX + s X 2 + (l+v) X ]=0 (71) 

and therefore 

sd s x +s X 2 +(l+v) X + 8 (t) = 0, (72) 

g(r) being an unknown function of t. For <f>(s, t), this implies: 

sd ss< f> + (1 + v)d s 4> + g( T )(f> = 0. (73) 

We now proceed by setting s = [h(T)y] 2 , with Ii{t) an arbitrary 
function of t, and also (f>{s) = [h{T)y\~ v t/f(y). In terms of i/f(y), 
Eq. d73l ) takes the following form: 

y 2 d yy ifr + yd y ifr + [Ag{T)h\ T )y 2 - v 2 ] ifr = 0. (74) 

If we now choose h(f) = — J=, we obtain: 

y 2 d yy <p + yd y i// + (y 2 - v 2 ) if, = 0, (75) 

which we recognize as the Bessel equation. The relevant so- 
lutions are xjj(y) = J v (y), and therefore 

0(s) = ^J v (2ViO>) (76) 

Alike in the case of the soft age, we again match the asymp- 
totic behavior of d69b for \s\ — > oo with the Green's function 
(l68l l at |z| — > 0. In this way (see the Appendix) we get that 

*(T)=£. 

C. Equation for the Cauchy transform of the orthogonal 
polynomials 

Following our previous paper [2], we conclude by briefly 
describing the evolution of the average of the inverse of the 
characteristic polynomial. In the case of our specific initial 
conditions (see after Eq. (TTOt). it is proportional to the follow- 
ing Cauchy transform of the time-dependent monic orthogo- 
nal polynomial: 

^ (Z ' T)= 2^Jo dX — g • 

(77) 

A straightforward (albeit tedious) calculation reveals that, 
given d55l >, the object defined by 

/^fcr) S Q V exp(^)/^(z,T) (78) 



satisfies: 

d T p V N (z, f) = -~ [zd zz p V N (z, T) + (1 + V) d z p V N (z, T)] , (79) 

a partial differential equation of exactly the same form as 
Eq. ( f55t . The difference resides in the associated initial con- 
ditions which, for the Cauchy transforms are singular at the 
origin. This results in a different choice of the corresponding 
solutions. In the case of the soft edge we recover the GUE 
result [2], i.e. Airy functions of the second kind. In the case 
of the hard edge, one gets Bessel functions of the third kind 
(Hankel functions), in agreement with Ref. [30] for the re- 
spective universality classes. 



V. CONCLUSIONS 

In this paper, we have studied a freely diffusing matrix 
of the Wishart ensemble. The Smoluchowski-Fokker-Planck 
equation for the stochastic evolution of its eigenvalues al- 
lowed us to construct, after a proper time rescaling, two par- 
tial differential equations. The first equation, which has the 
form of a generalized Burgers equation, describes the behav- 
ior of the associated Green's function in the limit of a matrix 
of infinite size. We have solved this equation with the method 
of complex characteristics for the particular case where ini- 
tially the matrix vanishes. We related the singularities of the 
characterstics to the occurrence of pre-shocks in the Burgers 
equation. Additionally, we translated this results into those 
describing chiral and anti-Wishart matrices and made con- 
nections to the method of free random variables. The sec- 
ond equation concerns the averaged characteristic polynomial 
which, for the chosen initial conditions, is equal to a certain 
time-dependent monic orthogonal polynomial of the Laguerre 
type. This partial differential equation, is an exact equation, 
for any matrix size, and it reduces to the generalized Burgers 
equation that has been obtained in the limit of infinite sizes. 
This equation lends itself to a scaling analysis in the vicin- 
ity of the edges of the large spectrum. We then recovered 
the asymptotic scaling functions describing the universal pre- 
shocs at the "soft edge", of the Airy type and, for the first time 
within the present approach, for the "hard-edge" singularity - 
of the Bessel type. This provides further evidence of the effi- 
ciency of the method. In a forthcoming publication |[l"8tl . we 
shall present a generalization of the derivation of the equation 
for the average characteristic polynomial that does not rely on 
the use of orthogonal polynomials. This will allow us to ex- 
plore other types of initial conditions leading to a new type 
of singularity. In particular, we derive for Wishart and chiral 
complex ensembles a new type of cusp singularity (general- 
ized Bessoid), which supersedes the known Pearcey cusp j27tl 
at the closure of the gap. 
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APPENDIX 

We start from matching the Airy edge. For |z| — > oo and 
< arg(z) < n, in the leading order j29tl 



Ai'jz) 
Ai{z) 



From ( f60b we get 



lim ri /R (z L /R(T) + N- 2/3 s, r) * A(t) + lim AT 1/3 

= A(T)+lim N^'h; 113 ^- 
x-xx Ai(x) 



(80) 



d s Ai(x(s)) 
Ai(x(s)) 



(81) 



where for simplicity we denoted the first term of the r.h.s. of 
d60l l as A(t). This expression has to match (for any argument 
of s except n) the corresponding limit for the resolvent, i.e. 



\l(z - z,)(z - z„) 
lim G(z) = A(t)- lim — — -, (82) 

z^z LIR 2rrz 



with z = Z m + N 2 ^ 3 s. An elementary calculation shows that 
Vfe - z L )(-s) _ 



and 



2rrz, 



y/(z R -z L )s 
2rrz„ 



: 2 r 3/2 T ' 



(83) 



(84) 



for the two edges, respectively, so that the above equalities 
impose the condition g(r) = 0. 

For the matching of the Bessel edge, let us note 12911 that 



(85) 



JAx) 

for |x| — > oo and < arg(x) < n. For the Bessel edge we have 

<I>{s,t) = s- vI2 J v {2^)s). (86) 

Therefore, for s — > oo and < arg(s) < 2n, a similar calcula- 
tion as the one above yields, asymptotically, 



fy{N- 2 s,T) = — -iN 



2t 



From the resolvent, we get 



lim G(z) = ^--i-L = ^-- iN-rL 

2=w- 2 ,n0 + 2r yrz 2t -\/ts 



(87) 



(88) 



which gives g(j) 



l 
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